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Abstract 

Dynamics of flexible ferromagnetic filaments in an external magnetic field is considered. We 
report the existence of a buckling instability of the ferromagnetic filament at the magnetic field 
reversion, which leads to the formation of a metastable loop. Its relaxation through three dimen- 
sional transformation of the configurations is observed experimentally and confirmed by numerical 
simulations. Bending modulus of the flexible ferromagnetic filaments synthesized by linking mi- 
cron size core-shell ferromagnetic particles with DNA fragments is estimated by comparison of the 
parameters of the loops observed in the experiment with theoretical calculations. Formation of 
the loop and its relaxation are characterized by the numerically calculated writhe number. The 
relaxation time of the loop allows us to estimate the hydrodynamic drag of the filament. 
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I. INTRODUCTION 



Ferromagnetic filaments are used by magnetotactic bacteria for the purpose of navigation 
in the magnetic field of the Earth jl, 2]. They may be created artificially by linking com- 
mercially available functionalized ferromagnetic microparticles with biotinized fragments of 
DNA [3]. Different phenomena are known for ferromagnetic filaments. Due to their flexibility 
ferromagnetic filaments orientate perpendicularly to the AC magnetic field if its frequency 
is high enough |3j. At magnetic field inversion they form loops [4|, which are metastable 
due to their migration to one of the ends of the filament. 

At the present moment the attention of researchers is focused on the creation of micro 
devices which mimic self-propelling microorganisms jlj]. It has been demonstrated that the 
chains of superparamagnetic particles may be driven by an AC magnetic field UQ. An es- 
sential feature of these developments is the symmetry breaking of the filament which causes 



its self-propulsion, 
of the filament 



"he symmetry breaking is achieved by linking some cargo to one end 
9], considering the defects in the chain of the magnetic particles j3, Q 
and inducing spatial symmetry breaking by the buckling instabilities of the magnetic rod 
[3, sj]. It is important to note that anisotropy of the hydro dynamic drag is necessary for 
the self-propulsion of the magnetically actuated swimmer 6MlO| . Neglect of anisotropy of 



the hydrodynamic drag invalidates different schemes pro 



Dosed for the magnetically actu- 



ated microswimmer (see for example the discussion in Due to their property of loop 



formation ferromagnetic filaments present an opportunity to consider a new mechanism of 
the actuation of the magnetic microswimmer 12J. By periodically inducing the buckling 
instability of the filament it is possible to create a microdevice which mimics the swimming 
of the biflagelate algae Chlamydomonos reinhardi 13( with the characteristic power and 



recovery strokes. Recently these algae have attained interest from the point of view of the 
study of synchronization phenomena [14|, the transition to multicellularity 15] and other 
interesting phenomena. 



Here the numerical algorithm based on 16] is developed for the numerical study of the 



stability of the loop in the 3D case. By numerical simulations it is shown that besides 
the instability due to the migration of loop there is a three dimensional instability. This 
conclusion is supported by experimental observation of the loop dynamics at magnetic field 
inversion. The formation of the loop and its relaxation are characterized by the numerical 
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calculation of the writhe number Wr. It is found that the characteristic lifetime of the 
loop depends on the magnitude of the perturbation - the closer the loop is to the planar 
configuration the longer its lifetime. 



II. MODEL AND ALGORITHM OF NUMERICAL SIMULATION 

Here we consider the buckling instability of the ferromagnetic filament induced by the 
reversion of the magnetic field. The physics of the process is the following. By applying a 
constant magnetic field a filament with the permanent magnetic moment orients along the 
field as a compass needle. At the magnetic field reversion the filament becomes unstable and 
relaxes to the energetically favorable state with the magnetization along the field. There 
are two pathways to how this may happen. One possibility is that both ends of the filament 
rotate in the same direction, for example clockwise. This so called S-like mode corresponds 
to the reorientation of the solid compass needle along the direction of the reversed magnetic 
field, which due to the flexibility of the filament occurs through intermediate S-like shapes. 
This mode of the filament relaxation does not have a threshold value of the magnetoelastic 
number Cm = MHL 2 /A, where M is the magnetization per unit length of the filament, 
2L is its length, and A is the bending modulus. Besides this the filament at Cm > it 2 /A 



171 ] has the relaxation mode when the both free ends rotate in opposite directions as is 



illustrated in Fig. 1. In the nonlinear stage development of this mode leads to the formation 



of the loop. It is metastable since as shown in the frame of the two-dimensional model J4] 
it relaxes by migration to one of the free ends of the filament. The characteristic relaxation 
time grows exponentially with the magnetoelastic number Cm and thus according to the 
two-dimensional model the loop is a long living metastable state. 

We illustrate here that there is a more effective process of loop relaxation through the 
third dimension. Its physics is the following, - for the loop the magnetization of its free 
ends is along the magnetic field and the bent fragment of the filament has a magnetization 
component antiparallel to the field. The magnetization of this fragment relaxes to the 
direction of the magnetic field by bending out of the plane of the loop. Necessary symmetry 
breaking for out of plane deformations arises due to random perturbations. 

Since by the formation of the loop the symmetry in the plane of the bending is broken 
then by inducing the periodic buckling of the filament in the AC magnetic field it is possible 
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to create a self-propelling micro device as it is illustrated in [121 ] . For the self-propulsion of 
the ferromagnetic filament the process of the loop relaxation through the third dimension 
may impose definite conditions on the time dependence of the magnetic field since the 
symmetry of the filament will be restored by the relaxation through the third dimension 
and the self-propelling motion will stop. We believe that this may be overcome by choosing 
a corresponding protocol for the time dependence of the applied magnetic field. 

The theoretical model for the dynamics of the ferromagnetic filament is developed in 



19| and is based on the expressions for the force F and torque T in the cross-section of the 
filament. It is assumed that magnetization is along the tangent vector of the filament (see 
Fig. 1). If the directions of the tangent do not coincide the torques arise and the filament 



may buckle. In the Frenet frame the relations for F and T read 

1 



T = -A-b , (2) 

where t is the tangent vector to the center line of the filament, n, b are the normal and 
binormal of the center line respectively, I is the arclength of the center line, R is the radius 
of the curvature of the center line and r its torsion. A characterizes tension arising due to 
the inextensibility of the filament, — Mt is the magnetization per unit length of the filament, 
H is the magnetic field strength. The contribution of the twist is neglected in (JTJ and (J2J). 
Different effects of the twist on the ferromagnetic filaments are considered in |19| . 

In the simplest case the dynamics of the filament is described in the Rouse approximation 

In the equilibrium F = if the external forces are absent. For the planar configuration 
introducing angle which the tangent makes with the field direction t = (cos ($), sin ($)), 
the condition F ■ n = leads to the equation 

A^ + MHsm($) = , (4) 

which in the case of the undamped ends at I = ±L is solved at boundary conditions 
dfi/dl \i=±l= 0. The solution reads 



•& = — 2arcsin (sin ( ~7^J sn y — ^ — > srn ) ) ) ' (^) 



Here ±i? m are the tangent angles at the free ends of the filament and sn(x, k 2 ) is the Jacobi 
elliptic function. The angle d m is found from the solution of the equation 

"3 

U r 



K 



(sin 2 (^))=V^, (6) 



here K(k 2 ) is elliptic integral of the first kind. Solution (jSj) for the maximal curvature of 
the filament at its center gives 

(4) =2sinf%)v^. (7) 



\RJ max V 2 

The dependence of the maximal curvature on Cm determined by the relation (j7|) is shown 
in Fig. 2. The Fig. 2 shows that at Cm > 10 the maximal curvature is well described by 



the asymptotic relation at dm — > tt (L/R) max = 2a/ Cm. Dependence shown in Fig. 2 may 
be used to determine the elastic properties of the ferromagnetic filaments by measurements 
of their curvature. The elastic properties of superparamagnetic filaments linked by PAA in 



the similar way are determined in 20]. 



Using the Frenet equations dt/dl = —n/R and dn/dl = t/R + rb the relation (CO) may be 
rewritten as follows 

F = -A^ + MH + At, (8) 

where A = A - Mt ■ H - 3A/2R 2 . Tildes further are omitted. 

The set of equations is solved numerically, discretizing the filament by a series of n + 1 
segments with the length h (a = ((h)^ 1 is mobility of segment) 

iT i = o l (f i + fr);(i = 0,...,n+l) , (9) 

here 

fi = -hAd 4 f/dl 4 | i5 fl = hd{Ai)/dl \i (i = 1, ...n - 1) ; (10) 
/ = -Ad 3 f/dl 3 | +MH; f' = At\ ; 
f n = Ad 3 f/dl 3 \ n -ME- f n = -Af U 

are the forces on the segments due to their bending, magnetic torque and tension. 

The stress At is found using the condition of the inextensibility tdv/dl = 0. It in the 
integral form gives 

= - / At^dl = At | v - At\ n v n + f v^pdl . (11) 



In the discrete form the condition (fill) reads 



71-1 



» + /X + ^/^ = 0. (12) 

i=i 

The condition of inextensibility is imposed by n constraints gt = {fk+i — r^) 2 = h 2 (k = 
0, ...,n — 1; j = 0, ...n). The n x cf(n + 1) Jacobian matrix with elements J = dgk/df*j is 
introduced. The relation f[T2"j) put in the matrix form as /' • v = is satisfied for allowed 
motions J ■ v = by /' = J T 7 (upper script T denotes transposed matrix, 7 is n-dimensional 
colon vector) [if]]. Equation of motion v = «(/ + /') accounting for constraints gives 
Jf + Jf = a~ x Jv = and thus 7 = —(J ■ J T )~ x Jf. As a result the equation of motion 
reads 

v = aPf, (13) 

where P = I — J T (J ■ J T )~ 1 J is the projection operator on the space of allowed motions 
(P 2 = P). 

Approximating d^f/dl 4 \f, (i = 1, n — 1) by the central differences, d 3 f/dl 3 |o= (r-j+rQ — 
2r*i)//i 3 and d 3 f/dl 3 \ n = — (r„_2 + r n — 2r n _i)//i 3 the stiffness matrix C, which determines 
the elastic stresses (/ = Cr + f m ) is constructed. The implicit scheme for the time step r 
then gives 

(J - arP ■ C)(f t+T - f t ) = arP ■ {Cf t + f m ) , (14) 

here /™ = MH; f™ = —MH; f™ = (i = 1, n — 1). The equation for r t+T — r t is easily 
resolved by inverting the matrix / — arP ■ C . 

At the end of each time stepthe filament is reshaped to satisfy the constraints according 
to the algorithm described in [16j. To put the equations in dimensionless form the following 
scales are introduced: time - ((2L) 4 /A - the characteristic elastic relaxation time, length - 
2L, and the elastic force A/(2L) 3 . The only parameter which controls the dynamics of the 
filament in this scaling is the magnetoelastic number Cm = MHL 2 /A. 



III. DYNAMICS OF FILAMENTS 



The synthesis of the ferromagnetic filaments by linking functionalized with streptavidin 
core-shell ferromagnetic particles with size d = 4.26 \im (Spherotech) with 1 kb long 
DNA fragments is carried out as described before {3J. Magnetic properties of the diluted 
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ensemble of these particles are measured by the vibrating sample magnetometer and are 
shown in Fig. 3. For the range of the magnetic field H < 200 Oe used in our experiments 
the magnetization curve in Fig. 3 gives the value of the magnetic moment of the particle 
m = 1.4 ■ 10~ 10 emu. The linear density of the magnetization M = m/d then is M = 
3.3 ■ 10 -7 emu. 

The protocol of the experiment is the following. The magnetic field is created by two coils 
in the space between two microscope slides at a distance of 50 fim. The coils are supplied 
by power supply- amplifier Kepco BOP 20-10M. A constant magnetic field is applied and the 
filament is oriented along the field lines. Then by switching the direction of the current in 
the coils in time less than 4 ms determined by the inductance of the coils the magnetic field 
is reversed. The dynamics of the configurations of the filament is registered by the video 
camera Mikrotron MC1363 at 50 frames per second. 

A sequence of several configurations of the ferromagnetic filament under the field inver- 
sion is shown in Fig. 4. It illustrates how the loop formed at an intermediate stage relaxes 
through transformation to a three dimensional configuration. The observed dynamics is 
compared with results of numerical simulations. Coordinates of the filament projection on 
x, y plane are obtained from the sequence of video frames by ImageJ. The contour length of 
the projection is calculated and the coordinates of the filament are smoothed by the poly- 
nomial of degree 11. Dimensionless coordinates of the filament are obtained by normalizing 
with the contour length of the initial configuration. The comparison of the results of the 
numerical calculations at with the experimental observations is shown in Fig. 5. To obtain 
the quantitative agreement of the parameters of the loop observed in experiment with nu- 
merical simulation it is carried out at Cm = 5.5. As initial configuration for the numerical 
simulation the initial configuration of the filament for the symmetry breaking inclined at 
angle a (tana = 10~ 4 ) to x, y plane is taken. Numerical simulation in agreement with the 
experiment shows that the filament makes a loop and does not relax to the straight configu- 
ration along the reversed direction of the field by the S-like deformation mode, which might 
develop at corresponding initial conditions The qualitative and for last 3-4 configu- 

rations quantitative agreement of the configurations obtained in experiment and numerical 
simulation is observed. The qualitative differences with the experiment start from the last 
configuration. In the numerical experiment filament through the third dimension rather 
fast relaxes, as it is illustrated in Fig. 6, to the straight configuration along the direction of 
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magnetic field. In the experiment the delay of the loop relaxation is observed and the end 
of the loop slides along its other leg for the time period 0.3 s, twice as large as the time 
necessary of the loop formation. We presume that it is due to the magnetic interaction of 
the end of the loop with the filament which is not taken into account in our model. The 
transition to the delayed relaxation stage of the loop is illustrated by the time dependence of 
the curvature of the filament in Fig. 7. The curvature is calculated by differentiation of the 
smoothing polynomial of degree 11. The error bars are estimated from 5-10 percent error 
in the coordinates of the filament due to their smoothing by the polynomial. To take into 
account the third dimension the curvature of the projection is corrected by the multiplier 
L p /L Q , where L p is the loop projection length and L is the length of the filament. Transition 
to the delayed relaxation stage takes place at 0.2 s and corresponds to the last configuration 
shown in Fig. 5. Relation (7) at Cm = 5.5 for the parameter (L/R) max gives the value 
4.3 not far from the value 4.85 corresponding to the final stage of loop formation shown 
in Fig. 5. For comparison with the theoretical curve of the curvature of the filament in its 
center given in a two dimensional model by the relation (J2J) in Fig. 8 we give the maximal 
curvature of the three dimensional filament and its curvature in the center obtained by nu- 
merical calculations at tana = 10 -3 . Corresponding curvature is given for the time moment 
when it is the maximal in a time sequence of the configurations for the given magnetoelastic 
number. We see that the curvature of the two dimensional model is between the maximal 
curvature and the curvature in the center of the filament. 

The value of the magnetoelastic number since the the linear magnetization density is 
known allows us to estimate the bending modulus of the filament A = 1.5 ■ 10~ 10 erg ■ cm. 
The characteristic time of the loop formation shown in Fig. 4 is 0.2 s and allows us to 
estimate the hydrodynamic drag coefficient £ = Airi] = 0.54 g ■ cm' 1 ■ s _1 which corresponds 
to the viscosity 0.043 g ■ cm' 1 ■ s' 1 approximately 4 times larger as the viscosity of water. 
This indicates that the hydrodynamic interaction of the filament with walls of the layer 
not taken into account in our model is important. Another indication about this are given 
by numerical simulations taking into account the anisotropy of hydrodynamic drag. It is 
incorporated introducing in the left side of equation (3) the matrix B' 1 = 7+(C|i — Cl)CZ 
Numerical calculation shows that at (<ji — CjJCJ 1 = — 1> which corresponds to the filament 
in unconfined conditions, relaxation to the straight configuration is more fast as it should 
be for the given value of the magnetoelastic number. The dynamics at — CjJCl 1 = 1 is 
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close to shown in Fig. 5. 

Besides the issues of hydro dynamic and magnetic interactions which have not taken into 
account by our model we should mention that the filament may possess different defects 
which may influence its dynamics |8|]. This may explain the quantitative differences of the 
filament evolution with the numerical simulation which may be seen in Fig. 5. 

The dynamics during the loop formation and its relaxation may be characterized by the 



writhe Wr 



2l| . which may be calculated according to the Fuller formula 22] 



2tt J l + cos(i?) ' v ' 

where •& is the angle which the tangent to the filament makes with the z axis. Time depen- 
dence of Wr for different inclination angles a of the plane of initial bending to the field at 
Cm = 6 is shown in Fig. 9. We see that the smaller the inclination angle the longer the loop 
lifetime. The sign of the Wr is determined by the sign of the inclination angle. Inclining 
the plane of initial bending in opposite direction Wr is positive. 

The loop formation and its relaxation at larger values of the magnetoelastic number has 
quantitative differences. As it is illustrated in Fig. 10 for Cm = 25 the loop formed at 
an intermediate stage has smaller radius, what is expected according to relation (J7|). We 
draw attention that the filament in Fig. 6 similarly to what is observed in the experiment 
(Fig. 3(3)) has configuration close to the ring what is not the case at Cm = 25 shown in 
Fig. 10, where the free ends are much longer. This confirms our estimate of the magnetoe- 
lastic number for the loop observed in the experiment. 



IV. CONCLUSIONS AND DISCUSSION 



To conclude it is shown that the loops of ferromagnetic filaments formed at magnetic 
field inversion are metastable and relax through the transformation to three dimensional 
configurations. Dynamics of the loop formation and relaxation is characterized by writhe 
calculations which show that the lifetime of the loop depends on the value of the initial per- 
turbation. Bending modulus of ferromagnetic filaments synthesized by linking the core-shell 
micron size ferromagnetic particles with the DNA fragments is estimated by the comparison 
of the characteristics of the loops observed in the experiment with the theoretical calcula- 
tions. 
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FIG. 1: Cartoon of a ferromagnetic filament at its bending in an external field. 
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FIG. 4: Loop formation and relaxation due to the escape into the third dimension at the inversion 
of the magnetic field. Time - 0(1);0.1 s(2); 0.16 s(3);0.28 s (4). H = 180 Oe,the length of the 
filament 2L = 74.6 /im. Field makes angle 7r/4 with the horizontal direction. 
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FIG. 5: Comparison of experimentally observed loop formation with numerical simulation (dashed 
lines). Dimensionless time starting from upper left corner increases with step 0.002 starting from 
0.002. Cm = 5.5. Observation time - increases by 0.02 s starting from 0.02 s. H = 180 Oe,the 
length of the filament 2L = 74.6 fim. Field makes angle 7r/4 with the x axis direction. 
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FIG. 7: Curvature of the projection of the filament at the loop formation and delayed stage of its 
relaxation before the escape into the third dimension. H = 180 Oe, the length of the filament 
2L = 88.5 fim. 
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FIG. 9: Time dependence of the Wr for several values of the inclination angle a of the plane of 
initial bending (tan (a) = 1CT 5 (solid line);10 -6 (long dashed line);10 -7 (short dashed line);10 -8 
(dotted line)). Lifetime of loop with smaller a is longer. Cm = 6. 
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FIG. 10: Dynamics of the loop relaxation as given by numerical simulation. Time in dimensionless 
units - 0.0034(1);0.0036(2);0.0038(3);0.0040(4). Cm = 25,tan(a) = 1(T 4 . Initial configuration 
y = e sin (nx/l)(l = 1 - (vre/2) 2 ;e = lO" 2 ). 
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